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ABSTRACT 



Aims. XTE J1701-407 is a newly discovered X-ray transient source. In this work we investigate its flux variability and 
study the intermediate long and short bursts discovered by Swift on July 17, and 27, 2008, respectively. 
Methods. So far, only one intermediate long burst, with a duration of «18 minutes and ten days later a short burst, 
have been recorded from XTE J1701-407. We analyzed the public available data from Swift and RXTE, and compared 
the observed properties of the intermediate long burst with theoretical ignition condition and light curves to investigate 
the possible nuclear burning processes. 

Results. The intermediate long burst may have exhibited a photospheric radius expansion, allowing us to derive the 
source distance at 6.2 kpc assuming the empirically derived Eddington luminosity for pure helium. The intermediate 
long burst decay was best fit by using two exponential functions with e-folding times of n = 40 ±3 s and T2 = 221 ±9 s. 
The bursts occurred at a persistent luminosity of L pcr = 8.3 x 10 36 erg s" 1 (~2.2% of the Eddington luminosity). For 
the intermediate long burst the mass accretion rate per unit area onto the NS was m « 4 x 10 3 g cm" 2 s~ 
energy released was -Ewst ~ 3.5 x 10 40 erg. This corresponds to an ignition column depth of 2/ign ^ U 
for a pure helium burning. We find that the energetics of this burst can be modeled in different ways, as (i) pure helium 
ignition, as the result of either pure helium accretion or depletion of hydrogen by steady burning during accumulation, 
or (ii) as ignition of a thick layer of hydrogen-rich material in a source with low metallicity. However, comparison of 
the burst duration with model light curves suggests that hydrogen burning plays a role during the burst, and therefore 
this source is a low accretion rate burster with a low metallicity in the accreted material. 
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Type I X-ray bursts are among the most evident signatures 
that testify the presence of a neutron star (NS) in low 
mass X-ray binaries. These bursts are thermonuclear 
explosions that occur on the surface of accreting NSs 
and are triggered b y unstable hydrogen and/or helium 
burning (see, e.g.. ILewin. van Paradiis fc Taaml . Il993t 
IStrohmaver fc BildstenI . 120061 for re views). Type I X-ray 
bursts were predicted theoretically bv lHansen fc Van Hornl 
(|l975r i. and sever al thousand burs t s have been observed 
to date (see, e.g.. ICornelisse et all l2003t iGallowav et all 
2006^ IChelovekov. Grebenev fc Sunvaevl l2006h . From the 
duration o f the bursts measured b y their decay parameter 
r (see e.g.. lGallowav et all l2006al) three main branches are 
distinguished: normal burs ts , intermediat e long bursts , and 
superbursts (see Fig. 7 in F alanga et all 120081 and refer- 
ences therein). These bursts can be described by different 
fuel types, accretion rates, and the y, therefore, also displays 
different recurrence times (e.g., [Strohmaycr & Bil dstenI . 
l2006tlCumming fc Macbethl 120041: [Galloway et alll2006af) . 
Thanks to the long Galactic plane scan carried out with 
BeppoSAX, INTEGRAL, and Swift a number of rare inter- 



mediate long bursts and superbursts have been observed. 
In most cases the rise of the burst is of w 1 s, whereas 
the decay is approximately exponential, with a duration 
of a few seconds for normal bursts, tens of minutes for 
intermediate long b ursts, and up to several hours for 
Kuulkersl 

" boost 



super bu rsts (e.g 



2004t iMolkov et al 



up 

2004; lin 't Zand et al 
in 't Zand et all 12005 



Chenevez etail . l2006l l2007t iFalanga et all . I2008D . The 



recurrence time of type I X-ray busts ranges from few 
hours to years, dep ending on the nuclear react io ns in- 
volved (see, e.g.. | Lewin. van Paradijs fc Taaml Il993t 
IStrohmaver fc Bildsteri . 20061 . for reviews). 
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In this paper, we report on an intermediate long burst 
from the X-ray transient XTE J1701-407. This source was 
discovered by the Rossi X-ray Timing Explorer (RXTE) 
during a routine Galact ic bulge scan on June 8, 2008 
()Markwardt et all l2008al) . At the time of the discovery, 
the source spectrum was best fit with an absorbed power- 
law model with photon index w 2.2 and a neutral ab- 
sorption column of « 3.4 x 10 22 cm" 2 , the measured flux 
was 1.4 x 10~ 10 erg c m" 2 s" 1 in the 2-10 keV band 
()Markwardt et all l2008afh Swift/XT'E follow-up observa- 
tions on June 11, 2008 found both spectral and source 
flux measurements consistent with the RXTE/PCA re- 
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Fig. 1. Flux history of the X-ray transient XTE J1701-407. 
Open squares represent fluxes derived from RXTE bulge 
observations in the period from March 31, to Septembre 2, 
2008. Fluxes derived from RXTE/PCA and Swift/XRT tar- 
get of opportunity and follow-up observations are marked 
with filled squares and stars, respectively. The arrow indi- 
cates the time of the type I X-ray burst. 



suits (|Deeenaar fc Wiinandsl . 12008ft . On July 17, 2008, the 
Swift/RAT camera detected a short flare consistent with 
the position of XTE J1701-407, an d ~ 97 s later Swi f t/XRT 
measured a decaying X-ray flux (|Barthelmv et all [2008). 
Based on this BAT data the X-ray spectrum was consis- 
tent with those observed for thermonuclear X-ray bursts, 
theref ore, this flare was associat ed to a most likely type I 
burst (|Markwardt et all l2008bD . Moreover, based on the 
Swift/XRT time reso lved spectrum during the flux decay, 
iLinares et al.l (|2008aD measured the cooling tail and con- 
firmed the flare to be an intermediate long type I X-ray 
burst. This led to the classifi cation of XTE J1701 -407 as 
a NS low-mass X-ray binary (|Linares et all 2008a, b). Ten 
days later, i.e., on July 27, BAT observ ed in the 15-150 
keV band a short bursts lasting « 10 s ()Sakamoto et all 
l2008h . Kilohertz quasi- periodic oscillations were observed 
at » 1150 Hz (see IStrohmaver. Markwardt fc Swankl . 
2008, for details). The most accurate position of the 
source was provided by Swift at aj2ooo=17 h 01 m 44!3 and 
fti2noo=-40°51 / 29"9 with an estimated accuracy of l'/7 
(jStarling fe Evanj . 12008ft . No infr ared counterpart candi- 
date was found at this position (|Kaplan fc Chakrabartvl 
12008ft . 



2. Data Analysis and results 

For the present study we used public available data 
from Swift and RXTE observatories. The dataset include 
RXTE target of opportunity observations, performed af- 
ter the discovery of the source (observation ID. 93444), 
as well as Swift follow-up observations (ID. 0031221001, 
0031221002, and 00317205001). The Swift data set includes 
also the intermediate long burst discovered on July 17, (ID. 
00317205000) and the short burst from July 27, 2008 (ID. 
00318166000). The total effective exposure time is 75.4 ks 
(30 pointings) and 1.1 ks for RXTE/PCA and Swift re- 
spectively. In Table Q] we report the detailed observations 
log. Data reduction of the RXTE Proport ional Counter 
Array (PCA; 2-60 keV, Uahoda et all . [1996) and the High 
Energy X-ray Timing E xperiment (HEXTE; 15-250 keV, 
iRothschild et all 11998ft was performed according to the 



default selection criteria for background subtraction, light 
curve, and spectrum extraction. We used only PCU2 data. 
We used Swift/XRT data in both Windowed Timing (WT) 
and Photon Counting (PC) mode. However, in the following 
we report only on WT data, since we found that PC data 
were strongly affected by pile-up and quasi-simultaneous 
observations carried out in WT mode were available. We 
reduced all the XRT data by using the version of the 
xrtpipeline (version 0.11.6) included in the Heasoft pack- 
age 6.4, and the latest calibration files available. Source 
events in WT mode were extracted from rectangle re- 
gions with widths of 40 pixels and heights of 20 pixelifl. 
Ancillary response files were generated with xrtmkarf and 
accounted for different extraction regions, vignetting, and 
point-spread function. For the Swift/RAT data analysis we 
used the batgrbproduct tool included in the Heasoft pack- 
age 6.4. Time resolved spectra and lightcurves (see Sec. l2.2p 
were extracted by using the batbinevt, batupdatephakw, 
batphasyserr, and batdrmgen tools. 

2.1, Persistent emission 

Spectra obtained from RXTE/PCA and Swift/XRT obser- 
vations were fit separately by using a photoelectrically- 
absorbed power-law model. The best fit parameters for 
the absorption column, Nn, are between 1.7 ± 1 and 
4.4 ± 1.2 x 10 22 cm -2 and a power-law index T = 2.0 ± 
0.1 — 2.7±0.3 the corresponding unabsorbed fluxes are be- 
tween l.l±0.07x 10" 10 and 8.3±0.6x 10~ 10 in the 2-20 keV 
energy band. In Table [1] we report all the persistent spec- 
tral parameters. All uncertainties are at a 90% confidence 
level. PCA spectra were extracted in the 2-20 keV energy 
band; the source was only detected at low significance in 
the HEXTE (15-50 keV band). A fit to the joined broad- 
band PCA/HEXTE (2-50 keV) spectrum did not lead to a 
significant improvement in the determination of the model 
parameters; therefore, we report in the following only on 
the PCA data. Swift/XRT spectra were extracted in the 2- 
10 keV band. All the measured unabsorbed fluxes were ex- 
trapolated to the 0.3-100 keV band by generating dummy 
responses (XSPEC version 11.3.2ag), and are shown in 
Fig. Q] We included in this figure also fluxes derived fr om 
RXTE bulge observation^ (ISwank fc Markwardtlf2001ft . In 
these cases, the conversion between RXTE count rate and 
flux was ob tained using the spectra r esults and the values 
reported bv lMarkwardt et all (|2008a[ ). 

The light curve and spectral analysis shows that 
XTE J 170 1-407 underwent flaring activity lasting a few 
days, during which the X-ray flux has varied by over a 
factor of 2-3. This variability can be ascribed to changes 
in the mass accretion rate or in the position of the 
source along its colour - colou r diagram (CCD, see e.g., 
lHasinger fc van der KlisJ . fl989f ). The latter possibility is in- 
vestigated in Fig. [2j The CCD was realized by using back- 
ground subtracted RXTE/PCA light curves with a 516 s 
time resolution. The soft colour is defined as the logarithm 
of the ratio between count rates in the energy bands 2.1— 
3.7 keV and 3.7-5.7 keV, whereas for the hard colour the 
energy bands 5.7-9.8 keV and 9.8-18.9 keV are used. No 



See 



also 



the 



XRT analysis manual 



at 



http://swift.gsfc.nasa.gov/docs/swift/analysis/xrt_swguide_vl_2.pdf 
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M. Falanga, et al.: Intermediate long bursts from XTE J1701-407 



3 



obvious transition between a hard to a soft state is ob- 
served. In order to reduce the errors on the colours, we 
generated a hard intensity diagram (HID) based on the net 
count rates in the 2.1-5.7 keV and 5.7-18.9 keV with a 516 
s time bin. Howevere, the source behaviour is not much bet- 
ter traced in the HID, because the statistical uncertainties 
are not significantly reduced along the two axes. In Fig. [2] 
we show the HID. Future observations are thus required in 
order to investigate further the origin of flux variations in 
XTE J1701-407. 
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Fig. 2. HID of XTE J1701-407. The hardness is the ratio 
of the count rates in the 5.7-18.9 keV to 2.1-5.7 keV and 
the intensity is in the 2.1-18.9 keV count rate. Each point 
corresponds to 516 s of integration time. 



2.2. Intermediate long burst light curves and spectra 

In Fig. El we show the Swift /BAT 15-150 keV (upper panel) 
and XRT 0.3-10 keV (lower panel) light curve of the inter- 
mediate long burst. The burst start time, is the time at 
which the BAT X-ray intensity of the source rose to 10% of 
the peak above the persistent intensity level. The XRT light 
curve starts with a delay of 133 s after the beginning of the 
burst as seen by BAT. The BAT light curve shows a slow 
rise tini£0 of w 45 s. The XRT decay time from the burst is 
best fit by using two exponential functions, with e-folding 
tim es of Tt = 44 ± 8 s and t 2 = 271 ±21 s, respectively (see 
also lLinares et al.l . l2008aD . The total duration of the burst, 
i.e. the time spent to go from and back to the persistent 
state, was of « 86 s and « 18 minutes in the BAT (15-150 
keV) and the XRT (0.3-10 keV) light curves, respectively. 

Type I X-ray bursts are produced by unstable burn- 
ing of accreted matter on the surface of the NS. The 
emission can be well described by a black-body radia- 
tion with temperatures, kT\^, in the energy range of a 
few keV. The energy dependent decay time of these burts 
is attributed to the cooling of the NS photosphere re- 
sulting in a gradual softning of the burst spectr u m. Fo r 
a review, see, e.g. iLewin. van Paradiis fc Taaml (|l993lh 
IStrohmaver fc Bildstenl (|2006l ).The time-resolved spectral 



The rise time is defined as the time spent between the start 
of the burst and the point at which the 90% of the peak burst 
intensity is reached. 
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Fig. 3. The intermediate long type I X-ray bursts detected 
from XTE J1701-407 on July 17, 2008. The time T ex- 
pressed in UTC corresponds to 17 h 31 m 55 s . The Swift /BAT 
(15-150 keV) and XRT (0.3-10 keV) light curve are shown 
with a time bin of 5 s. The data gap of f» 133 s in the XRT 
light curve is due to the time elapsed between the BAT 
trigger and the follow-up observation. The dashed line in- 
dicates the background level measured « 3000 s after the 
burst end. 



analysis of the burst was carried out by using BAT and XRT 
data in the 15-35 keV and 0.3-10 keV bands, respectively. 
Burst spectra were well fit by a simple photoelectrically- 
absorbed black-body, BB, model. For all these spectral fit 
the Nh value was frozen at 3.1 x 10 22 cm™ 2 derived from the 
pre-burst persistent emission (see Table [1]). The inferred bb 
temperature, fclhb, apparent bb radius at 6.2 kpc (see Sec. 
13. 1|) . Rhb, and bolometric luminosity are shown in Fig. [5] 
and in Tabled The burst fluence is calculated from bolo- 
metric fluxes, Fboi! these correspond to the observed 2- 
10 keV Swift/XTR fluxes extrapolated to the 0.1-100 keV 
energy range. The peak flux, -F pe ak, is derived from the BAT 
15-35 keV light curve spectra and extrapolated to the 0.1- 
100 keV energy range. All the measured unabsorbed fluxes 
were extrapolated to the 0.1 100 keV band by generating 
dummy responses (XSPEC version 11.3.2ag). This is well 
justified for the XRT data since the black-body tempera- 
ture is well inside the spectra bandpass. However, during 
the main peak a black-body te mperature reaches a max- 
imun between 2-3 KeV (see e.g.. iFalanga et all l2008i and 
references therein), and therefore extrapolating the black- 
body spectra outside the BAT bandpass. For the BAT (15- 
35 keV) burst peak black-body spectral best fit we found 
fcT bb = 2.68 ± 0.21 (x 2 = 0.74, 7 d.o.f.). We fixed an up- 
per and lower boundary black-body temperature of 2 and 3 
keV, respectively, living the normalization as a free parame- 
ters. Note the lower boundary is also consistent to be higher 
as the first XRT black-body temperature of fcTbb ~ 1-8 (see 
Table [3]). We found an unacceptable fit by using a 2 keV 



black-body burst peak temperature (% 2 



5.1, 8 d.o.f.), 



whereas a slightly better fit is obtained by using a temper- 
ature of 3 keV (x 2 cd = 1-22, 8 d.o.f.). In Fig. g] we show 
the BAT data with the different black-body models from 
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Table 1. Log of RXTE and Swift observations and best fit spectral parameters of the persistent emission. An absorbed 
power-law model is used to fit these spectra. 



Instrument 


Date 


Exp. JVh 


r 


FluX2_iokcV 


FluX2_20keV 


Fl ux 0.3-100keV Xred 




MJD 


ks 10 22 cm -2 




erg cm -2 s _1 


erg cm -2 s _1 


erg cm -2 s _1 



Swift b 


54628. 


.66385 


0. 


.4 


4 0+°' 7 


2 


9+0.3 
• z -0.3 


2. 


8 ± 0.3e-10 


3 


.7 


±0. 


5e-10 


10.2 ± 2e-10 





.8 


RXTE 


54633. 


.24637 


2. 





3.4+ 2 2 


2 


4+0.1 

•^-0.2 


4. 


.5 ± 0.4e-10 


5 


.6 


±0. 


,4e-10 


17.5 ± 7e-10 


1 


.2 


RXTE 


54634. 


.42452 


2. 


5 




2 


4+0.1 

•^-0.2 


3. 


.4 ± 0.2e-10 


4 


.3 


±0 


,2e-10 


12.8±2e-10 





.7 


RXTE 


54635. 


.27545 


2. 





q q+2.1 
°- J — 2.2 


2 


9+0.1 
• z -0.2 


2. 


7±0.2e-10 


3 


.0 


±0 


,2e-10 


9.4 ± 1.6e-10 





.8 


RXTE 


54638. 


.80878 


3. 


5 


o 4+O.8 


2 


9+0.1 
z -0.1 


1. 


.9 ±0. le-10 


2 


.6 


±0. 


.le-10 


6.8 ± 1. le-10 





.9 


RXTE 


54639. 


.85563 


3. 





o 4+0.7 


2 


q+0.1 
°-0.1 


2. 


.1 ±0.1e-10 


2 


.8 


±0. 


.le-10 


7.8 ± 1.4e-10 


1 


.2 


RXTE 


54640. 


.90230 


3. 


.5 


o r+0.7 
°-°-0.7 


2 


9+0.1 
z -0.1 


1. 


.9 ±0. le-10 


2 


.0 


±0 


.le-10 


7.1 ± 1. le-10 


1 


.2 


Swift b ' c 


54664. 


.62684 


0. 


6 


q 1+0.4 

J - i -0.4 


1 


q+0.2 

a -0.2 


4. 


.9 ± 0.3e-10 


7 


.2 


±0 


,6e-10 


18.3 ± 3e-10 


1 


.1 


RXTE 


54665. 


.10546 


3. 





3.7+ 2 2 


2 


fi +0.2 
°-0.1 


4. 


.9 ± 0.6e-10 


5 


.9 


±0. 


,7e-10 


20.0 ± 4e-10 


1 


.6 


RXTE 


54669. 


.24527 


1. 


,8 


2.9+ 2 2 


2 


E+0.2 
J -0.2 


4. 


.8 ± 0.8e-10 


5 


.8 


±0 


,9e-10 


18.0±3e-10 


1 


.6 


RXTE 


54670. 


.99462 


3. 


3 


3.9+ 2 2 


2 


c+0.1 
d -0.2 


5. 


5 ± 0.7e-10 


6 


.5 


±0 


,6e-10 


23.2 ± 4e-10 


1 


.5 


RXTE 


54672. 


.97638 


1. 


7 


4.1± 2 ' 2 


2 


7+0.2 
• ' -0.2 


4. 


.9 ± 0.5e-10 


5. 


.8 


±0 


,6e-10 


21.5±5e-10 


1 


.3 


RXTE 


54674. 


.73749 


1. 


.8 


4.2+ 2 2 


2 


7+0.2 
• 1 -0.3 


4. 


.6 ± 0.7e-10 


5 


.5 


±0 


,6e-10 


20.8 ± 5e-10 


1 


.6 


Swift b ' d 


54674. 


.94054 


0. 


1 




2 


u -0.4 


7. 


.3 ± 0.5e-10 


8. 


.3 


±0 


,6e-10 


23.3 ± 4e-10 


1 


.4 


RXTE 


54677. 


.47471 


3. 


5 




2 


fi +0.3 
D -0.2 


4. 


5 ± 0.7e-10 


5 


.3 


±0 


,8e-10 


18.3±3e-10 


1 


.2 


RXTE 


54678. 


.59119 


3. 


3 


4.4+ 2 - 2 


2 


q+0.3 
°-0.1 


4. 


.3 ± 0.7e-10 


5. 


.0 


±0. 


,8e-10 


17.5 ± 4e-10 


1 


.5 


RXTE 


54679. 


.63915 


3. 


2 


4.ll 2 ' 2 


2 


7+0.2 
• 1 -0.3 


4. 


0±0.5e-10 


5 


.1 


±0 


,6e-10 


14.9 ± 5e-10 


1 


.3 


RXTE 


54680. 


.29712 


2. 


1 


2.911' 2 


2 


c+0.2 
d -0.3 


3. 


7±0.5e-10 


4 


.3 


±0 


,6e-10 


13.2 ± 3e-10 


1 


.5 


RXTE 


54681. 


.86101 


3. 


5 


3.6+i 2 


2 


q+0.2 
°-0.2 


3. 


.3 ± 0.5e-10 


3 


.7 


±0 


,6e-10 


14.7±3e-10 


1 


.1 


RXTE 


54682. 


.71212 


7. 





o q+0.9 
z - J -0.9 


2 


c+0.1 
J -0.1 


1. 


5 ± 0.04e-10 


1 


.9 


±0 


,05e-10 


6.5 ± 1.4e-10 


1 


.1 


RXTE 


54683. 


.89045 


3. 


.5 


9 9+I.2 
2 


2 


9+0.1 
• z -0.2 


0. 


.95 ± 0.07e-10 


1 


.2 


±0 


,5e-10 


3.4 ± 0.6e-10 





.8 


RXTE 


54684. 


.93786 


1. 


8 


22 +o.9 


2 


9+0.2 
• z -0.2 


0. 


.89 ± 0.03e-10 


1 


.4 


±0. 


,3e-10 


3.2±0.5e-10 





.9 


RXTE 


54685. 


.80434 


1. 


6 


9 n+1.0 
z - u -0.6 


2 


9+0.1 
z -0.1 


0. 


.84 ± 0.06e-10 


1 


.1 


±0 


,07e-10 


3.0±0.6e-10 


1 


.2 


RXTE 


54686. 


.31211 


2. 


.5 


i q+0.9 
r.O_ .4 


2 


q + 0.1 
•°-0.1 


0. 


.94 ± 0.06e-10 


1 


.2 


±0 


,6e-10 


3.5±0.8e-10 





.8 


RXTE 


54687. 


.04156 


0. 


7 


4.0+ 2 


2 


c+0.2 
d -0.2 


1. 


.2 ±0. le-10 


1 


.5 


±0. 


,2e-10 


5.2 ± 1.2e-10 





.7 


RXTE 


54687. 


.35989 


2. 


7 


9 7+1.2 
^•'-1.2 


2 


q+0.1 
°-0.1 


1. 


1 ±0.06e-10 


1 


.5 


±0 


,06e-10 


4.4 ± 0.9e-10 


1 


.3 


RXTE 


54687. 


.43638 


2. 


2 


9 7 +1.3 
z -'-1.3 


2 


q + 0.1 
•°-0.1 


1. 


1 ±0.07e-10 


1 


.4 


±0 


.08e-10 
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Unabsorbed flux. " Swift/XRT was operating in WT mode. c ' Intermediate long and short pre-burst persistent emission, 

respectively. 



the different fits. In this case the best fit values and the 
extrapolated peak flux are acceptable values in the order of 
20% wich is inside our BAT extrapolation error box. 

In Table [5] we report all the measured bust parameters. 



Table 2. Parameters of the Intermediate long burst. 



^eak (erg cm" 2 s" 1 ) 8.2 ± 2 x lO" 8 
fl (erg cm" 2 ) 7.6 ± 0.3 x 10" 6 

r ee /b/F peak (sec) 92 ± 22 
7 ee F p c cr5 /f peak 2.2 ±0.5 x 1Q- 2 



"Unabsorbed flux (0.1-100 keV). 
c Unabsorbed persistent flux F pors = 
s" 1 (0.1-100 keV). 



b Fluence (0.1-100 keV). 
(1.8 ±0.2) x 10~ 9 erg cm -2 



From the time-resolved spectral results we converted 
the light curve count rates to flux and found also in this 



case, that a double exponential function are required to fit 
the intermediate long burst decay. The derived e-folding 
time are t\ w 40 ± 3 s and T2 ~ 221 ± 9 s, respectively 
with a Xicd = 1-1 (f° r 163 d.o.f.). In this fit we included 
also the first two BAT data points. Note, using only the 
XRT data points we found the same results within the er- 
ror bars. A single exponential function is found to be in- 
adequate with Xr C d = 5-2 (166 d.o.f.). In order to compare 
the decay tail with the intermediat e long X-ray burst from 
2S 0918-549 a nd SL X 1737-282 (|in 't Zand et all 120051 : 
iFalanga et all 120081). and with t he decay cooling model 
from lCumming &: Macbe"th1 (|2004t ) we fitted the data with 
a power-law, and found an index of T = —2.14 ± 0.02 with 
a x^ ed = 1.7 (167 d.o.f.). In Fig. [S]we show the double ex- 
ponential best fit (upper panel) and in a log-log scale the 
power-law best fit (lower panel) to the data. The double 
exponential fit is statistically preferred over an power-law 
fit by the F-test probability of 1.2 x 10~ 14 . Similar values 
of ti = 44 ± 7 s and T2 = 218 ± 55 s or power-law in- 
dex r = —1.97 ± 0.09 can be obtained from the fit to the 
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Fig. 4. Swift/B AT (15-35 keV) data of the intermedial long 
peak black-body spectra and best fit model, kThh ~ 2.7, 
(solid line). Upper limit, kThh — 3 keV, acceptable fit 
(dashed-line) and lower, kThh = 2, keV unacceptable fit 
(dotted-line). 
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Fig. 5. Evolution of the spectral parameters, as inferred 
from Swift/BAT (first 86 s) and Swift /XRT (833 s) obser- 
vations. The bolometric luminosity is calculated by assum- 
ing a distance of 6.2 kpc, see Sec. 13.11 The bottom panel 
shows the xt 



/2 cd values for these fits. 



bolometric luminosity reported in Fig. [5J Also in this case 
the double exponential fit is statistically preferred over an 
power-law fit by the F-test probability of 6.6 x 10 -3 . 

The decay time determined with the light curve in unit 
of count rates, is restricted to the XRT 0.3-10 keV energy 
band. This decay time does not take into account the differ- 
ent rise and cooling effect for different energy bands during 
the outburst, where converting the rates to the bolometric 
flux or luminosity using the time-resolved spectral results 



a more realistic decay time can be obtained. Therefore, in 
the following we consider only the T\ ~ 40 s and T2 rj 221 
s or the power-law index —2.14. In Table [2] we report the 
t = /b/Fp Ca k = 92 ±22 s which is only valid if one exponen- 
tial function describe the decay tail. In our case this value 
is consistent within the error bars if we consider the double 
exponential decay as; /b/F P oak = e(Ti — T2) + t~2 = 71 ± 12, 
where e = F x /F pea ^ = 0.83 and F x is the peak flux for n. 
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Fig. 6. Top panel: We show the BAT and XRT intermediate 
long burst light curve with the double exponential best fit 
curve. Bottom panel: The same data presented in a log-log 
scale with the power-law best fit model. 



Table 3. Fit parameters of the time-averaged BAT and 
XRT burst spectra. 
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2.3. Short burst light curves and spectra 
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Time (sec) 



Fig. 7. Top panel: The short type I X-ray burst detected 
from XTE J1701-407 on July 27, 2008. For comparison of 
the rise time and duration we show also the intermediate 
long burst (see also Fig.EJ). The Swift/BAT (15-150 keV) 
light curves are shown with a time bin of 2 s. 

In Fig.[7]we show the Swift/BAT 15-150 keV short burst 
light curve. The time T expressed in UTC corresponds 
to 22 h 31 m 19 s . The rise time is 1 ± 0.5 s where the total 
duration is « 8 s. To determine the rise time and burst 
duration we rebined the light curve to 1 s. For this short 
burst the XRT light curve starts with a delay of 180 s after 
the beginning of the burst as seen by BAT. In the XRT 
data (0.3-10 keV) the burst has not been detected, so the 
total duration of this burst has to be considered < 180 s. 
In Fig. [7] we show the short burst, and for comparison we 
plot also the intermediate long burst seen by BAT. 

Given the short duration and statistic of the burst light 
curve we were not able to study a time-resolved spec- 
tral analysis. Here we report the burst spectrum analy- 
sis integrated over 6 s in order to measure as accurate 
as possible the peak flux. Burst spectra were well fit by 
a simple black-body model. The inferred bb temperature, 
kTbb = 3.9 ± 0.5 keV with an unabsorbed bolometric peak 
fluxes, Fboi = 3.8 ± 1.1 x 10" 8 (0.1-100 keV). The burst oc- 
curred at a flux persistent level of 2.3 ±0.4 x 10 -9 erg cm~ 2 
s _1 . This short burst occurred at a comparable persistent 
emission level as the intermediate long burst. 



3. Discussion 

3.1. Source distance, persistent flux, and accretion rate 

When a burst undergoes a photosphcric radius expan- 
sion (PRE), the source distance can be determined based 
on the assumption that the bolometric peak luminos- 
ity is saturated at the Eddington limit, ^Edd , (e.g . , 
iLewin. van Paradiis & Taaml 119931 : iKuulkers et "all . [2003) . 
During the PRE episode, while the bolometric luminosity 
remains constant at the Eddington value, the high energy 
flux may display a double peak profile and/or a delay in 
the rise time (e.g.. IKuulkers et all 120021 : iGallowav et all 



l2006al iFalanga et all l2007t ). TheBAT light curve shows a 
slow rise time of « 45 s, which is typically observed at high 



energy in intermediat e long helium bursts with PRE (e.g. , 
Kuulkers et"aT1 . l2002t iMolkov et all . 120051 : IFalanga et all 



20081) . 



Since Swift/XRT missed the first 133 second of the burst 
from XTE J1701-407 we cannot directly determine whether 
this burst underwent a PRE. This issue cannot be resolved 
with BAT time-resolved spectral analysis due to limited 
statistics (see Fig. [5]). However, the comparison with other 
PRE bursts showing this slow rise time at high energy leads 
us to take the observed profile in the BAT light curve as an 
evidence for a PRE during the first « 50 s of the interme- 
diate long burst, in which case this should correspond to 
the timescale for the photosphere to fall back to the neu- 
tron star surface. For the short burst no conclusion can be 
drawn as to whether a PRE event occurred (see Fig. [7]). 

Another possibility is that the burst has an intrinsi- 
cally long rise time. For ex ample, the mixed H/H e bursts 
observed from GS 1826-24 ([Galloway et all l2004h do not 
show PRE, but h ave a rise lasting f or 10 seconds set by 
hydrogen burning (jHeger et alll2007t) . The duration of the 
bursts from GS 1826-24 is much less than the long burst 
from XTE J1701-407. It is not clear whether a long duration 
mixed H/He burst could have a 50 second rise (see discus- 
sion of light curves in section 3.5). The fact that helium 
burning can be extremely rapid, whereas hydrogen burning 
involves slow weak interactions means that the rise time 
is longer when hyd rogen is present , and PRE much less 
likely. For example, iFuiimoto et al.l (|1987f ) derived a criti- 
cal helium fraction necessary to achieve PRE. 

Assuming a bolometric peak luminosity equal to the 
Eddington value for a He type I X-ray burst (i^Edd ~ 
3.8 x 10 38 ergs -1 , as empirically derived bv lKuulkers et all 
2003), we obtain the source distance of d = 6.2lJ'g kpc. 
For comparison, the theoretical valu e of this distance (e.g., 
ILewin. van Paradiis fc Taaml . I1993D found by assuming a 
He atmosphere and canonical NS parameters (1.4 solar 
mass and radius of 10 km), is 5.5^0 g kpc. Note that the 
source could be closer if the peak luminosity of the burst 
was lower than the pure helium Eddington limit. For exam- 
ple, we cannot rule out that the burst did not show PRE, 
in which case the peak luminosity could have been sub- 
Eddington. Alternatively, taking the peak luminosity to be 
the Eddington luminosity for solar composition (Xq = 0.7) 
gives a distance of ~ 4 kpc. In the following we consider 
d 6.2 kpc as a fiducial distance, and comment on how 
our conclusions would change if the source was closer. 

The best fit to the 2-20 keV persistent emission spec- 
trum of XTE J 170 1-407 required an absorbed simple power- 
law model with T « 2.1. Assuming a distance of 6.2 kpc, 
the estimated intermediate long burst pre-burst persistent 
unabsorbed flux between 0.1-100 keV, F x 



pcrs 



1.9 x 10" 



erg cm s , translates into a bolometric luminosity of 



L pors « 8.3 x 10 36 erg s 



or w 2.2%L Ed d- 
The local accretion rate per unit area is then given by 



4TrR 2 m(GM/R)/(l 
4.0 x 10 3 g cm" 2 s" 1 



R 



11.2 km 



M 



1.4 M, 



6.2 kpc 



1 + z 
1.26 



1) 



A convenient unit for accretion rate is the Eddington ac- 
cretion rate. Here, we define the local Eddington accretion 
rate as mEdd = 1-8 x 10 5 g cm -2 s _1 , the local accretion 
rate onto a M — 1.4 Mm. and R = 11.2 km neutron 
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star that gives an accretion luminosity equal to the em- 
piric ally derived Eddington luminosity 3.8 x 10 38 erg s" 1 
from lKuulkers et all (|2003f h This gives m/mEdd = 2.2%. 

3.2. The energy, ignition depth, and recurrence time of the 
long burst 

The observed energy of the long burst allows us to estimate 
the ignition depth. The measured fluence of the burst is 
fb = 7.6 x 10~ 6 erg cm" 2 , giving a burst energy release 
E buist = 4ird 2 f b = 3.5 x 10 40 ergs (d/6.2 kpc) 2 . The ignition 
depth is given by £ b urst = 47ri? 2 y ign Q Iluc /(l + z), or 



?/ig n = 1.8 x 10 9 g cm 2 



Qn 



6.2 kpc 
-l 



R 



l + z 

1.6 MeV/nucleon J V 11 - 2 km J V I- 26 



(2) 



The value of Q nuc ~ 1-6 MeV corresponds to the nuclear 
energy release per nucleon for complete burning of helium 
to iron group elements. Including hydrogen with mass- 
weighted mean mass fraction (A) give s Qnuc ~ 1.6 + 4 (A) 
MeV/nucleon (| Galloway et all l2004f ). wh ere we include 
losses due to neutrino emission following iFuiimoto et alj 
(1987). For (X) = 0.7, the solar composition value, Q nuc = 
4.4 MeV/nucleon, and yi gn — 6.5 x 10 s g cm' 2 . 

At an accretion rate of 4 x 10 3 g cm" 2 s — 1 , the re- 
currence time corresponding to a column depth of yi gn — 
1.8 x 10 9 g cm" 2 (pure helium composition) is At = 
(z/ign/?ti)(l + z) = 6.6 days, or for y- lgll = 6.5 x 10 8 g cm -2 
(solar composition) is At = 2.4 days. These derived recur- 
rence times are independent of the assumed distance. The 
intermediate long burst, as the first observed burst, occured 
40 days after the detection of XTE J1701-407. The effec- 
tive exposure time on the source from the beginning of the 
outburst to the intermediate long burst was 0.34 days, and 
0.14 days between both bursts. The elapsed time on the 
source is thus too short compared to the theoretically de- 
rived recurrence time to allow us to get an observational 
measurement of the recurrence time. 



3.3. Theoretical comparison with ignition models 

To try to understand the nuclear burning processes respon- 
sible for the long burst, we compare the observed properties 
with type I X-ray burst ignition mo dels. The ignition con- 
dition s are calculated as described in lCumming fc Bildstenl 
(|2000l ). but we take R = 11.2 km, M = 1.4 M Q , and the en- 
ergy release in hot CNO burning^ Eh = 6.0 x 10 18 erg g _1 . 
The results are shown in Table 2] for several different choices 
of accreted hydrogen fraction Xq, accretion rate rh and flux 
from the crust Qb- The ignition conditions are calculated by 
modeling the temperature profile of the accumulating fuel 
layer, and adjusting the thickness of the layer until a condi- 
tion for thermal runaway is met at the base. These models 
do not include the effects of previous b ursts on the igni- 
tion conditions (e.g. IWooslev et aTll2004f) . and we have not 
included gravitatio nal sedimentation which is important at 
low accretion rates (jPeng et all l2007f ) . In the following, we 



4 He re we adopt the value of Eh from Wal lace fc Wooslevl 
lll98lD. which includes neutrino losses. These were not included 
bv lCumming fc Bildstenl IpOOOl ). 



consider three posibilities: 1) accretion of pure helium, 2) 
hydrogen-rich matter at solar, and 3) low metallicity. 

3.3.1. Pure helium accretion 

For pure helium accretion, the accumulating fuel layer is 
heated from below by a flux rhQb from the crust, where Qb 
is the energy per accreted nucleon flowing outwards from 
the crust. We set the accretion rate at the observed value 
w/rh-Edd = 2.2%, and adjust the value of Qb until we find 
ignition at the inferred column for pure helium burning 
Uign = 1.8 x 10 9 g cm" 2 (and therefore obtain the observed 
burst energy). A flux from below Qb — 0.5 MeV/nucleon 
matches the observed burst energy at 2.2% Eddington. We 
also include models with Qb — 0.3 and 0.7 MeV/nucleon in 
Table [4] to show the sensitivity of the ignition depth to the 
amount of heat from below. Note that since the combina- 
tion rhQb sets the total flux heating the layer, an increase or 
decrease in Qb keeping rh fixed is equivalent to increasing 
or decreasing rh with Qb fixed. Therefore, if the inferred 
accretion rate is smaller by some factor, the value of Qb 
needed to match the observed energy will be larger by the 
same factor. Similarly, a source distance less than 6.2 kpc 
can be accommodated by increasing Qb- 

The outwards flux expected from the crust dep ends on 
accre tion rate and the core neutrino emissivity (|Brownl 
I2000D . ranging from rhQb ~ 0.1 MeV per nucleon at high 
accretion rates up to rhQb ~ 1.5 MeV per nucleon at low ac- 
cretion rates. The value depends on how much of the total « 
1.5 MeV per nucleon heat released in th e crust by pyenonu- 
clear and electron capture reactions (|Haensel fc Zdunikl . 
1990, 2008) is conducted into the core and released as neu- 
trinos compared to being conducted outwards. The num- 
ber of 0.5 MeV/nucleon that we infer here is quite reason- 
able for an accretion rat e of 2.2% Eddington; fo r example, 
the models calculated bv lCumming et ail(|2006l ) for persis- 
tently accreting sourc es have Qb = 0-3 to 0.9 MeV/n ucleon 
at this accretion rate. lGallowav fc Cumm ing ( 2006b) mod- 
elled the X-ray bursts observed from SAX J1808.4-3658 and 
found Qb ~ 0.3 MeV/nucleon for rh ss 3% mEdd- 



3.3.2. Accretion of hydrogen rich matter with solar 
metallicity 

Next, we consider that the source is accreting hydrogen 
rich matter with the solar hydrogen fraction Xq — 0.7 and 
a metallicity similar to solar, with mass fraction of CNO 
elements .Zcno = 0.02. In that case, the hydrogen burns at 
a fixed rate during accumulat ion of the fuel l a yer, b y the 
beta-limited hot CNO cycle of iHovle fc Fowled |l965D . The 
hydrogen depletes at a column depth 



lid 



1.5 x 10 8 g cm" 



0.022m Ed d 



^CNO 

0.02 



Ao 
0.7 



(3) 



smaller than the inferred ignition depth, so that a thick 
layer of pure helium accumulates and ignites beneath a 
steady hydrogen burning shell. Therefore, even though so- 
lar composition material is accreted, the mean hydrogen 
fraction at ignition in the fuel layer is very small, giving 
Qnuc close to the value 1.6 MeV/nucleon for pure helium. 

For a CNO mass fraction Zcno — 0.02, the observed 
burst energy is obtained for an accretion rate three times 
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Table 4. Type I X-ray burst ignition conditions 



Model 


(%m Edd ) 


ZcNO 




Q b 








Qnuc 




At 

(days) 


Pure helium accretion 1 * 


1 


2.2 


0.02 





0.5 


1.8 1.4 








1.6 


3.4 


6.5 


2 


2.2 


0.02 





0.3 


7.4 1.2 








1.6 


14 


27 


3 


2.2 


0.02 





0.7 


0.84 1.5 








1.6 


1.6 


3.1 








Accretion 


hydrogen rich material 










4 


2.2 


0.02 


0.7 


0.5 


0.13 2.0 


0.38 


0.07 


3.1 


0.49 


0.48 


5 


0.69 


0.02 


0.7 


0.5 


1.6 1.4 


0.01 





1.64 


3.2 


19 


6 


2.2 


0.001 


0.7 


0.1 


0.67 1.7 


0.62 


0.54 


4.1 


3.3 


2.5 


7 


2.2 


0.001 


0.7 


0.5 


0.53 1.8 


0.63 


0.43 


4.1 


2.7 


2.0 



a Models 1, 5, and 6 provide a good match to the observed burst energy of 3.5 x 10 ergs. In addition, models 1 and 6 have an 
accretion rate that matches the value inferred from the persistent luminosity. 

b We define iriEdd = 1-8 x 10 5 g cm" 2 s" 1 , the local accretion rate onto a 1.4 Mq, R = 11.2 km neutron sta r that gives an accretion 
luminosity equal to the empirically derived Eddington luminosity 3.8 x 10 38 erg s _1 from iKuulkers et al.l (|2003f ). 
c The hydrogen mass fractions are: in the accreted material Xo, at the base of the layer at ignition Xf,, and the mass- weighted 
mean value in the layer at ignition {X}. 

d Note that the ignition conditions for pure helium accretion do not depend on the choice of Zqno ■ 



lower than inferred from the observed X-ray luminosity 
rh/mEdd — 0.69% (model 5 in Tabled]). At an accretion 
rate of 2.2% Eddington, the burst energy is a factor of seven 
too small (model 4). Reducing the distance to the source 
does not help because it would change both the inferred m 
and burst energy by the same factor. 



3.3.3. Accretion of hydrogen rich matter with low metallicity 



In the third scenario, we consider that the material is hy- 
drogen rich, but with a low metallicity. Model 6 has a burst 
energy close to the observed energy at the inferred accre- 
tion rate of 2.2% Eddington, with a CNO mass fraction 
Z = 10~ 3 , roughly 10% of the solar metallicity. In this 
case, the amount of hot CNO burning is reduced substan- 
tially, so that hydrogen permeates the entire fuel layer at 
ignition. Hot CNO burning still operates at a low level and 
causes some preheating of the fuel layer. The hydrogen in- 
creases the amount of nuclear release during the burst, giv- 
ing Qnuc ~ 4 MeV/nucleon, and an ignition depth three 
times smaller than for the pure helium case. The low metal- 
licity ignition is less sensitive to Qb than the pure helium 
ignition. Model 7, which has the same conditions as model 
6, but with Qb — 0.5 rather than Qb = 0.1, has a burst 
energy within almost 20% of the observed value. Similar to 
the pure helium ignitions, a closer distance can be accom- 
modated by varying Qb- 

One caveat regarding the low metallicity ignition models 
is that there can be substantial heating of the accumulating 
fuel layer because of nuc lear reactions a s sociat ed with the 
ashes of previous bursts. IWooslev et al.l (|2004t ) found that 
the burst behavior at accretion rates « 0.1 Eddington was 
insensitive to the metallicity of the accreted material due 
to this effect. However, they found that at lower accretion 
rates there was go od agreement with th e ignition models 
presented here (see IWooslev et all . 12004 Table 9). 



3.3.4. Summary of ignition models 

We present three ignition models in Table [4] that fit the 
observations. Models 1 and 6 have the correct burst en- 
ergy and accretion rate. They correspond to accretion of 
pure helium (model 1), for which the layer is heated by the 
outwards flux from the neutron star crust, and for accre- 
tion of hydrogen rich material with low metallicity (model 
6) , for which a low level of hydrogen burning preheats the 
layer during accumulation, but the hydrogen fraction at ig- 
nition is substantial and makes a significant contribution to 
the burst energetics. Third, model 5 has the correct burst 
energy, but at an accretion rate three times lower than ob- 
served. Given the uncertainties in translating the observed 
X-ray luminosity to accretion rate, it seems worthwhile con- 
sidering this model further. In this model, the accreted com- 
position is hydrogen-rich with a solar metallicity. This leads 
to depletion of the hydrogen by the hot CNO cycle and the 
build-up of a thick layer of pure helium beneath the hydro- 
gen shell. 

There is a fourth possibility, which is that the source is 
accreting hydrogen-rich material, but the hydrogen burns 
unstably in a series of short flashes. The helium produced 
in the short flashes builds up and make a pure helium layer 
that ignites to give the long burst. This is similar to model 
5, but with unstable rather than stable hydrogen burning. 

One way to distinguish the various possible explana- 
tions for the burst energetics would be a recurrence time 
measurement: as shown in Table |4j the different scenarios 
predict differences in recurrence time. In fact, this is equiv- 
alent to a measurement of the a parameter^ for the bursts, 
which would indicate the fuel type (e.g. a ~ 40 for solar 
hydrogen abundance as found for example in GS 1826-24, 
iGallowav et~aTll2004l ). 

While this paper was in preparation, a similar analysis 
of the intermediate burst from XTE J170 1-407 was carrie d 
out and reported in the preprint by iLinares et alJ (l2008bTl . 
The bolometric peak flux for the long burst and therefore 

5 The quantity a is defined as the ratio of the total energy 
emitted in the persistent flux to that emitted in a burst, a = 
.FpersAt//b, where At is the time interval between two bursts. 
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the limits on the source distance are consistent with the 
values we obtain here. Ho wever, the persistent luminosity 
and burst energy given in iLinares et al.l (|2008bf ) are a fac- 
tor of two smaller than our values, being restricted to the 
energy range 1-50 keV, whereas we have included a bolo- 
metric correction. These values are used in the interpre- 
tation of the long burst, resulting in an accretion rate of 
m = 1.9 x 10 3 g cm~ 2 s _1 = 1.1% mEdd and a burst en- 
ergy of Eb = 1.6 x 10 40 ergs, both a factor of two smaller 
than the values we find in this paper. Repeating the igni- 
tion calculations presented earlier for this lower accretion 
rate, we find that accretion of solar composition material 
(Xq = 0.7) gives a burst energy of Eb = 1.5 x 10 40 ergs 
for Qb = 0.1 MeV/nucleon, and Eb = 1.0 x 10 40 ergs for 
Qb = 0.5 MeV/nucleon. Therefore, accretion of solar com- 
position naturally explains the burst energetics for these 
values of Eb and m. For pure helium accretion at this rate, 
we find that Qb = 1.5 MeV/nucleon is required to achieve 
a burst energy of 1.5 x 10 40 ergs (lower values of Qb result 
in a deeper ignition and more energetic burst). Whereas 
this is a larger value than expected at this accretion rate, 
it is within the range of the t otal energy released in the 
crust (jHaensel fc Zdunikl . [2008) . In addition, since the flux 
heating the layer is oc rhQb, the requirements on Qb can 
be relaxed if the true accretion rate is larger than assumed 
he re. Therefore, we find that the energetics argument given 
in ILinares et al.l (|2008bl) against explaining this burst as 
pure helium is overstated. 

3.4. Constraints from the light curve 

The shape and duration of the light curve offer another 
way to determine the composition of the fuel that burns 
during the burst. While we do not have models available 
with exactly the same burst energy , the low accretion rate 
models from IWooslev et al.l (|2004f ) are within a factor of 
two to four in energy and accretion rate, and so we compare 
these models with the observed light curve. These models 
have M = 3.5 x 10~ 10 M® yr , corresponding to a local 
accretion rate m = 1% mEdd- 

First, we compar e the observed l ight c urve with burst 
2 from model zm of IWooslev etal! (|2004f ) (solid curve m 
Fig.[8|). This model is for accretion of hydrogen rich (Xq w 
0.7) matter with a low metallicity, Z = 10~ 3 , and so is 
similar to model 6 in Tabled] The burst has a total energy 
release of 2.0 x 10 40 ergs, just less than a factor of two 
smaller than the observed burst. The recurrence time is 3.0 
days, and ignition column y- lgn = 3.6 x 10 8 g cm -2 . The 
ignition column is just less than a factor of two smaller 
than the observed burst. The peak luminosity of this burst 
is close to the Eddington luminosity for solar composition, 
suggesting that the distance may be closer than the 6.2 kpc 
assumed here. A closer distance would bring these light 
curves into better agreement. The model light curve has 
a steep decline at late times, somewhat steeper than the 
observed decline. An extra factor of two in ignition column 
would give the correct burst energy and lengthen the model 
burst light curve, bringing it into better agreement with the 
observed burst. Another point to note is that this burst has 
a slow rise time, lasting for several seconds, as exp ected for 
a low helium mass fraction (|Fuiimoto et alJ . H987| ). 

T he dashed curve in F ig. [8] shows burst 3 from model 
Zm of lWooslev et al.l (|2004l ). This model has hydrogen-rich 



matter with solar metallicity Z — 0.02, and so is similar 
to model 5 of Table [5] The hydrogen burns away in a thin 
shell, leaving a pure helium layer that ignites. This burst 
has a total energy release of 8.0 x 10 39 ergs, a factor of 
4 smaller than observed. The recurrence time is 4.5 days, 
and ignition column depth 5.6 x 10 s g cm~ 2 . The burst 
reaches close to the pure helium Eddington luminosity. The 
duration of the burst is shorter than the observed burst by 
a factor of w 5-10. The slope of the decay is shallower than 
the observed slope. The rise time of this burst is very fast, 
a fraction of a second, in contrast to the much slower rise 
of the hydrogen-rich burst (model zm; solid line in Fig. [8]) . 

We also include some cooling models calculated follow- 
ing [Cumming!^^acbitb| (|2004T ). For a given ignition col- 
umn, an energy release per gram of 1.6 MeV per nucleon is 
deposited in the layer, as would be appropriate if the he- 
lium burned to iron group elements at each depth at the 
start of the burst. We have also computed models for a 
lighter ash and corre spondingly smaller energy deposition 
(|Wooslev et al.l 120041 find that the burning does not go all 
the way to iron group in their model zM), but the differ- 
ences are small, and this does not change our conclusions. 
The cooling of the layer is then followed, with the flux from 
the surface limited to the Eddington luminosity for pure 
heliurr0. We show two examples in Fig. [5) The fi r st has 
y ign = 6 x 10 s g cm" 2 to match the IWooslev et al.l (2004) 
model Zm burst. The second has yi gn = 2x 10 9 g cm -2 as 
needed to get the observed burst energy (models 1 and 5 in 
Table 3]) . The shape of the cooling models agrees well with 
the model Zm light curve, and agrees within a factor of two 
in the cooling timescale. Even allowing for this factor of two 
in the y- lgn — 2 x 10 9 g cm -2 model, the cooling is faster 
than the observed light curve. 

In summary although the models from IWooslev et al.l 
(2004) are not at exactly the same ignition conditions as 
implied by the observations of XTE J1701-407, our com- 
parison suggests that pure helium ignition at the inferred 
ignition column depth (with or without a small overly- 
ing hydrogen burning shell) have cooling times that are 
shorter than the observed light curve. On the other hand, 
a hydrogen-rich composition throughout the layer, as ex- 
pected for low metallicity, produces a longer lasting light 
curve more consistent with the observed cooling time. This 
conclusion does not depend on the assumed distance to 
the source. Both helium and hydrogen-rich burst models 
reach the Eddington luminosity (either the pure helium or 
solar composition Eddington luminosity respectively), and 
therefore could explain the PRE suggested by the similarity 
between the observed BAT light curve and other interme- 
diate long bursts (see section 3.1). If the slow rise time is 
intrinsic to the burst and not due to PRE, the hydrogen 
rich model is preferred as the presence of hydrogen leads 
to a much slower rise time than for pure helium (Fig. [5]). 
The double exponential nature of the decay may also argue 
for hydrogen burn ing during the burst. The burst profiles 
from GS 1826-24 ([Galloway et al.l . l200l are well-fit by a 
double exponential decay. Further modeling is required to 



6 As noted by IWooslev et al.1 (|2004D . the shape of the light 
curve as the luminosity begins to drop below the Eddington 
luminosity is likely not accurately reproduced by these models, 
which do not follow the outer layers in detail. We will improve 
our treatment of this in future work. However, we expect that 
the late time cooling is not sensitive to this. 
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study the expected burst profiles of hydrogen-rich bursts 
produced as a result of low metallicity accretion. 
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Fig. 8. Model light curves compared with the observed 
light curve. The solid cu rve shows burst 2 from model zm 
of IWooslev et alJ (|2004h . This burst has an ignition col- 
umn and energy a factor of 2 smaller than the observed 
burst. The dashed cu rve shows burst 3 from model Zm of 
Woo slev et al.l ([2004D . This burst has an ignition column 
a factor of three to four times smaller than the observed 
burst. The observed bolometric flux has been converted 
to luminosity using a distance of 6.2 kpc. Redshift cor- 
rections have been applied to the theoretical light curves, 
with 1 + z = 1.26. The do tted curves show two coolin g 
models calculated following ICumming fc Macbeth! ((2004), 
for (left to right) y- lgn = 6 x 10 8 g cm -2 corresponding 
to Woosley et al. (2004) model Zm (dashed curve), and 
2/ign = 2 x 10 9 g cm~ 2 as needed to explain the observed 
burst energy (models 1 and 5 in Tabled]). 



3.5. Origin of the short burst 

We have focused on the long duration burst, which has a 
well-measured fluence and therefore energy. For the short 
burst, we can only place an upper limit on its fluence. At 
low accretion rates, unstable ignit ion of hydrogen can give 



rise to short duration burst s (e.g., Strohmayer fc Bildstenl . 



[2fM iChenevez etafl 12551 . ILinares et alT p008b) suggest 



that this is the origin of the short burst from XTE J1701- 
407, and that either (i) hydrogen-ignited short bursts are 
producing the helium fuel for the long burst, or that (ii) 
the source is accreting close to the boundary between un- 
stable and stable hydrogen burning, and stable hydrogen 
burning produces the helium for the long burst. The com- 
parison to ignition models and model light curves that we 
have made earlier suggest that hydrogen survives to the 
ignition depth, implying a low metallicity in the accreted 
layer. Unfortunately, this would presumably make a ther- 
mal instability driven by CNO burning less likely. 

For the low metallicity model or for pure helium ac- 
cretion, another explanation is that the local accretion rate 
was higher at the time of the short burst. The ignition depth 
is very sensitive to the base flux or e quivalently to the prod - 
uct rhQb (see for example Fig. 8 of lin 't Zand et aLll2005h . 
The persistent flux at the time of the short burst from 



XTE J1701-407 was slightly larger than at the time of the 
long burst, even at 30%, not enough to cause a significant 
reduction in the ignition column depth. We note that short 
bursts were observed from the intermediate long X-ray 
burster 2S 0918-549, a suspected ultracompact binary and 
therefore accreting hydrogen-deficient matter. Therefore, it 
is not clear to us that the observation of the short burst 
rules ou t pure helium accreti on in XTE J1701-407, as ar- 
gued bv ILinares et all (|2008bD . 

4. Conclusions 

We have compared the observed properties of the long dura- 
tion burst from XTE J1701-407 with models of type I X-ray 
burst ignition conditions and light curves. We showed that 
the observed burst energy could be understood as (i) pure 
helium ignition, either as a result of pure helium accretion 
or of depletion of hydrogen by steady burning during ac- 
cumulation, or (ii) as ignition of a thick layer of hydrogen- 
rich material with low metallicity. Comparing with model 
light curves, we find that the pure helium ignitions cool 
faster than observed. On the other hand, a hydrogen rich 
layer gives a longer duration light curve with a steep de- 
cline in the tail of the burst, better matching the observed 
light curve. Therefore we suggest that the intermediate long 
burst from XTE J1701-407 was powered by unstable burn- 
ing of a thick layer of hydrogen rich matter with low metal- 
licity. Long X-ray bursts caused by pure helium ignitions 
benea th a hydrogen shell have been i dentified, as for exam- 
ple bv lGallowav fc Cummind (|2006bf ) who argued that this 
was happening in SAX J1808. 4-3658. But to our knowledge 
XTE J1701-407 would be the first example of a source that 
shows long bursts driven by a thick layer of hydrogen-rich 
material. The bursts from GS 1826-24 are believed to be 
powered by rp-process hydrogen burning, giving long « 100 
second tails, but the ignition depth in those bursts is an or- 
der of magnitude smaller than inferred for the long burst 
from XTE J1701-407, so that hydrogen can survive until he- 
lium ignition, even for solar metallicity. At the low accretion 
rate in XTE J1701-407, this is not the case: low metallicity 
is required to reduce the rate of hot CNO burning and al- 
low hydrogen to survive until helium ignites. This implies 
either that this source is a burster accreting low metallicity 
H-rich material at a low rate, or another possibility is that 
heavy elements are able to sediment ou t from the accumu - 
lating layer at this low accretion rate (jPeng et all l2007t ). 
reducing its effective metallicity. If so, future studies of this 
source could be used to test the physics of sedimentation 
at low accretion rates. 
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